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The current-source density (CSD) analysis is a widely used method in brain electrophysiology, but this 
method rests on a series of assumptions, namely that the surrounding extracellular medium is resistive and 
' uniform, and in some versions of the theory, that the current sources are exclusively made by dipoles. Because 

^ ■ of these assumptions, this standard model does not correctly describe the contributions of monopolar sources 

| or of non-resistive aspects of the extracellular medium. We propose here a general framework to model electric 

• fields and potentials resulting from current source densities, without relying on the above assumptions. We 

, develop a mean-field formalism which is a generalization of the standard model, and which can directly incor- 

porate non-resistive (non-ohmic) properties of the extracellular medium, such as ionic diffusion effects. This 
formalism recovers the classic results of the standard model such as the CSD analysis, but in addition, we pro- 
vide expressions to generalize the CSD approach to situations with non-resistive media and arbitrarily complex 
multipolar configurations of current sources. We found that the power spectrum of the signal contains the sig- 
nature of the nature of current sources and extracellular medium, which provides a direct way to estimate those 
properties from experimental data, and in particular, estimate the possible contribution of electric monopoles. 

ir 



PACS numbers: 87.19.La, 87.17.Aa 



I. INTRODUCTION 



> 

On 

The current-source density (CSD) analysis [1, 3] is a method consisting of estimating the underlying current sources 
from a series of recordings of the extracellular electric potential. This method is widely used in neuroscience, and applies 
well to layered structures of the brain, such as cerebral cortex, hippocampus or cerebellum [1]. The CSD analysis is 
based on the "standard" model of electric potentials in biological tissue [1, 3, 4], which rests on the hypothesis that the 
extracellular medium is resistive (ohmic) and uniform. Other influences, such as ionic diffusion, are assumed to play a 
negligible role on the propagation of the electric field. 

•i-H . 

. Based on this set of hypotheses, the equation that determines the electric potential at macroscopic scales (~ 50 or 

^ \ more) is given by: 

v • (o- c vv) = o- L y 2 v = ^ , (i) 

at 

where <x e is the electric conductivity of the extracellular medium. This expression can be obtained by applying the 
differential law of charge conservation and Ohm's law. The term is interpreted as the volumic density /,„ of current 
sources. This equation forms the basis of the CSD analysis method [1-3, 5]. 

According to Eq 1, the electric potential V would only depend on electric conductivity and not at all on electric per- 
mittivity. However, Poisson's law in a homogeneous medium (eV 2 V = -p) implies that V will be twice smaller for twice 
larger s with the same charge distribution, so it is paradoxical that permittivity is not taken into account in CSD analysis. 
Moreover, according to Eq 1, the electric potential is determined solely by the charge conservation law, and independently 
of Poisson's law, which is contradictory with Gauss' law in Maxwell equations. 
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If we take Gauss' law into account, we can write 

ot s 

The general expression for p is: 

p(x, t) = p(x, 0)e _7 Mw 

where t M w = ■§; is the Maxwell-Wagner time of the medium. However, t M w is usually considered as negligible (typical 
values for biological tissue are t mw w 10~ 10 s, with cr = 0.3 S /m and s « 10~ 10 F/m) such that the current source density 
must be approximately zero, which is paradoxical. One way to resolve this paradox is to consider that the Maxwell- 
Wagner time is not negligible, or that electric parameters display strong spatial variations. However, such conditions 
contradict the hypothesis that the medium is resistive, and lay outside the domain of validity of Eq. 1 because in this case 
the impedance of the extracellular medium is complex in Fourier space [6, 7]. 

Despite this paradox, the standard model seems to apply relatively well to media such as brain tissue [1,3]. This model 
has, however, the drawback that it cannot be used to determine the validity of the hypotheses it is based on. Moreover, 
there is no clear definition of microscopic or macroscopic levels, and consequently, it is difficult to include possible 
frequency dependencies that could result from different physical phenomena at intermediate (mesoscopic) scales, such as 
ionic diffusion or membrane polarization [7]. 

In the present paper, we introduce a more general formalism which does not rest on the hypotheses of classic CSD 
where the medium is hypothesized to be uniform and resistive, which also supposes that the electric parameters are 
constant and independent of frequency. The goal of this new formalism is to provide generalized expressions for CSD 
analysis in non-resistive media. Our aim is also to provide a theory which is general enough to enable testing different 
hypotheses concerning the nature of current sources and the electrical properties of the surrounding extracellular medium, 
which could then be directly estimated from experimental data. 



II. GENERAL THEORY 



In this section, we derive a mean-field theory of the electric field and potential resulting from current densities in 
biological tissue, by staying as general as possible. 



A. Definitions and scales 



In the generalized formalism presented below, we will define the current sources from conductance variations. We will 
assume that the differential law of charge conservation holds in a given domain D, without defining a current source density 
per unit volume. We assume that the conductance variations in cellular membranes (especially around synapses), within 
domain D, are the principal origin of the extracellular electric field. This assumption is more realistic and biological 
compared to the "classic" assumption which is based on current source densities per unit volume because real current 
sources are caused by the opening or closing of membrane conductances in neurons. Note that these two different points 
of view can be complementary if we assume that the volumic density of current source I m is 




In this case, the two points of view are mathematically equivalent. This is the reason why we wrote the source term as 
-finEq.l. 

An important assumption of the present formalism is that all observable phenomena can be modeled by fields which 
are twice differentiable (class C 2 ). While most fields will obey this criterion, it will considerably simplify the mathemat- 
ical analysis to its simplest expression (commuting spatial and temporal first-order derivatives). In mean-field physics, 
by virtue of the Stone-Weierstrass theorem, it is always possible to make a uniform-convergence approximation of the 
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observable phenomena by a mean-field model of class C 2 . Indeed, because the mean-field of a discontinuous field of first 
kind is necessarily a continuous field (the primitive of a discontinuous function of first kind is continuous), this restriction 
to class C 2 fields will not limit the applications of the theory developed here. Moreover, the theory will remain general 
because most fields are class C 2 in practice. It will not apply to very particular models, such as fields involving surfaces 
with infinitely small thickness, current sources without volume, fields that necessitate infinite energies [22] 

Finally, the formalism developed below is only valid for well-defined ranges of spatial and temporal scales. We will 
consider scales greater than 1 pm (about 300 times the size of ions such as K + or Na + , including solvatation). This scale 
is chosen large enough for classic electromagnetism theory to apply without ambiguity (although in principle it can apply 
to scales down to a few nm). We will define as "microscopic", scales of the order of 1 pm, while "macroscopic" scales 
will be of the order of 50 /mi or more. We will also consider the typical range of frequencies of electrophysiological 
signals, up to 10 kHz, for which the quasistatic approximation is valid. 



B. Mean-field Maxwell theory 



We start from Maxwell equations where we consider the spatial averages of the fields and electric parameters, which 
will be denoted here by brackets < ... >. The spatial average is made over some reference volume, which is invariant. 
Because of the scale invariance of Maxwell equations (e.g., see [8], chapter 4), the spatial averages of electric field 
E, electric displacement D, magnetic induction 3 and magnetic field ft are linked by the following linear operatorial 
equations: 

v- < 6 > = < p fm > v- < 3 > = o 

(2) 

Vx<£> = Vx<//> = <)>+' 2 ^ 

where < j*> and < pf ree > are the spatial averages of the current density and free charge density, respectively. 

These equations allow one to find the general regularities that all models must satisfy. For example, the laws of energy 
conservation and momentum conservation can be deduced from this set of equations [8]. In particular, by using the 
relation V • (V x C) = (which is in general true for all vectorial fields of class C 2 ), one can deduce the differential law of 
charge conservation: 

d <D> _ dV- <D> _ d< p free > 

V • (Vx < A >) = V- < f > +V • — = V- < f > + = V- < f> + £ = (3) 

at ot at 



However, because the above equations relate the spatial averages of interaction fields (E,D,B,H) with the spatial 
averages of the two matter fields {f,pf ree ), it is necessary to complete them with a specific physical model to apply them to 
a given biological medium. This specific model must allow measuring spatial averages at a scale which is determined by 
the measurement method (type of electrode for example). Thus, the measurement system determines a minimal reference 
volume, which necessarily implies to use a mean-field formalism. 

In general, (for all media of class C 2 ), the fields <E>,<D>,<H> and < B > are linked by the following general 
equations: 

<D*>(r,t) = J ^ < e > (r,T,E,H) < E > (r,t - t) dr + < C > (r,t, E, H) 

(4) 

<B>{r,t) = £ + J <h>(?,t,E*,H) <ft>(f,t-T) dr 

where fi and e are respectively the absolute magnetic permeability and absolute electric permittivity tensors. Here, we 
have defined < D* > = < D > + < C > where < C > is the source field resulting from conductance variations. 
Note that in classic electromagnetism, one defines the electric displacement relative to vacuum permittivity by < 
z5 w > = So, < E a > + < Pa > (in frequency space) [23], which expresses the fact that the polarization field is 
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proportional to the electric field through electric susceptibility < Xoj > (< P<o > - < Xu >K E w >)■ It follows that 

< til > = < Soo >< £ a > + < P u > + < 6 a >= < So >< £ a > + < 6 a >. 

Considering Maxwell-Gauss' law (V- < D > = < pf ree >) and the definitions of the interaction fields imply the 
following relations between charge density and interaction fields: 

V<D*> = <p{ ree > 
V-<eoo><£> = < p{ ree > + <p A cond > + < p bo "" d > 

(5) 

V- < P > = -< p bo "" d > 



v- < C > = -<p h cond 



> 



where < p " > represents the average variation of free charge density produced by conductance variations, and 
< p{ ree > is the average free charge density which does not result from membrane conductance variations. Note that the 
divergence of the field < C > depends on the exact mechanism of conductance variation. If this mechanism does not 
produce monopoles, then this divergence is zero. Also note that the field < C > is generally assumed to be independent of 
the field < Z? >, which is a valid assumption for biological media in general, except if ephaptic interactions must be taken 
into account. 

We also have 

< f> f) = I < o* >(r,T,M,fl) <i>(f,t-r)dr+<D> V < p free > . (6) 

*J — oo 

where cr e is the electric conductivity and < D > is the mean ionic diffusion tensor [24]. We define < D > as follows: 

N N 

< d > v < P s r >=Yu <Di> w< P' ree > 

i=l i=\ 

where the sums run over the different ionic species. Thus, we can write that the part of current density caused by 
concentration changes equals < D > V < pf ree >[25]. In expression (6), we have separated the current produced by 
ionic diffusion from the current produced by other physical causes such as Ohm's effect, polarization, etc. Note that this 
separation was made here for simplicity, but it is also possible to integrate diffusion effects in the expression of the mean 
conductivity (see Eq. A4 in Section A 1). 

It is important to note that the first term in the righthand side of Eq. 6 is not exclusively due to Ohm's law (which 
relates to energy dissipation), but can reduce to it in some cases [7, 10]. In general, Eq. 6 gives a time-dependent 
electric conductivity (or frequency-dependent in frequency space), which is not the case for Ohm's law in general (see 
Appendix A 3). Also note that the integrals in Eqs. 4 and 6 can be seen as convolution products relative to time, in which 
case they take the form of a simple product in frequency space. 

Thus, according to this theory, it is sufficient to model the physical and geometrical nature of the extracellular medium 
by using electromagnetic parameters and diffusion coefficients to simulate the interaction fields when the current sources 
are known - this is usually called the forward problem. Inversely, we can also deduce the physical characteristics of the 
sources from the knowledge of the electromagnetic parameters and interaction fields, as well as their spatial and temporal 
variations - this is known as the inverse problem. 

Finally, the integrals in Eqs. 4 and 6 must satisfy the causality principle, according to which the future cannot determine 
the present state of the system. For example, the value of electric field < E > at time t + \At\ must not influence the value of 
electric displacement < D > at time t. Thus, the causality principle determines a supplementary constraint on the possible 
values of tensors < /j >, < s > and < <x p >, which limits the number of possible mathematical models of the extracellular 
medium. For instance, as detailed below in Section II F, this principle imposes mathematical relations between the electric 
parameters < s > and < cr e >, which are called Kramers-Kronig relations for linear media. 

The set of equations above define a mean-field formalism in which Maxwell equations are formulated with spatial 
averages. In the next sections, we consider different approximations to this formalism. 
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C. The quasi-static approximation in mean field 



5 



The first approximation to the Maxwell equations is the quasi-static approximation, which consists of de-coupling 
electric and magnetic variables. In general, the time variation of < 3 > produces an electric field < M > (Lenz-Faraday 
effect), the electric and magnetic variables are coupled in Maxwell equations (Eq. 2) by the following expression: 

d < B > 

Vx<£> = (7) 

at 



It was shown that for biological media and current sources, the Lenz-Faraday effect is negligible [11]. In such condi- 
tions, we can write: 



Vx < E > = (8) 



Under this quasi-static approximation, the electric and magnetic variables are de-coupled in Maxwell equations, and 
the electric field obeys: 

<e> = -v<y> 



This approximation is also called the a priori quasistatic approximation, by opposition to the a posteriori quasistatic 
approximation, which consists of finding the general solution of Maxwell equations and later de-couple the electric and 
magnetic variables (see details in [11]). Although these two approximations are not strictly equivalent, we will only 
consider the a priori quasistatic approximation in the remaining of this paper. 

According to this approximation, Maxwell equations simplify to the following expressions: 

.free 



V- < D > = < p' ree > 



Vx < E > 



(9) 



where the current density < j* > is linked to the electric field < E > by Eq. 6. Note that, contrary to the static cases 
(electrostatics and magnetostatics), the fields <£?>, <£)>, <j*> and < p$ ree > are here space and time-dependent. 



In the following, we consider the complex Fourier transform 



f 

<J —a 



X(t)e~'°"dt , X(t) 



2n J_ 



X,„ e ,L "dcj 



where a> — 2nf. Note that because of the linearity of the spatial average, we have < X > w =< X w >. 
Applying the complex Fourier transform to Eqs.4, 6 and 9 leads to: 



V • (< e w > V < W u >) 
{ V • (< ff-J > V < V a >) 



free 



> + V- < c,„ > 



ico < piZ" > - io) V- < 6 m > + V • (< D > V < pi ree >) 



(10) 



The first of these equations is Poisson's law in mean-field. Although in some cases (electrostatics), the Poisson equation 
is sufficient to determine the solution of the system, it is not sufficient in the quasistatic case, and the second equation 
is necessary to close the system. This second equation is the differential law of charge conservation (in the presence of 
diffusion) and takes into account the time variations of the electric potential 



Bedard & Destexhe, revised (September 15, 201 1) 6 
By multiplying the first equation of Eqs. 10 by id), and adding the result to the second equation, leads to: 

v 2 < v bJ > + v<Tm> • v < y„ >= — !— v ■ (< d > v < P s r >) (ii) 

where y w =< crj > +/w < s u > is the admittance of the extracellular medium. This equation is general and can be 
used to calculate the extracellular potential in an extracellular medium with arbitrarily complex properties (i.e., when the 
electric parameters depend on frequency and space). It is a generalization of expressions obtained previously [6, 7]. The 
righthand term accounts for ionic diffusion. 

Because Maxwell equations are scale invariant, the expression above (Eq. 11) is valid at all scales. Like in any mean- 
field approach, the spatial scale can be chosen according to the scale of the phenomenon that needs to be modeled, as well 
as the physical size and distance between electrodes. For example, in the case of CSD of mammalian cerebral cortex, one 
must consider scales of the order of 50 //m to resolve the field produced by each cortical layer. 

Note that in the quasistatic approximation, the explicit dependence of the electric field on magnetic permeability p 
completely disappears. However, there can still be an implicit dependence through H in nonlinear media, because the 
electric field does not depend explicitly on magnetic induction anymore. 

D. The quasistatic approximation at larger scales 

At small scales (» 1 pm), biological media such as the cerebral cortex are far from homogeneous and isotropic. The 
electric parameters can display large variations, for instance between fluids and membranes. However, at larger scales 
(« 50 pm), such media can be considered as homogeneous and isotropic. In such a case, the tensors < e al > and < crj > 
can reduce to scalar quantities. Note that the fact of considering larger scales suppresses the directional dependence of 
the propagation of currents by a statistical equivalent, without changing the frequency dependence produced by physical 
phenomena at small scales. The transition from small scales to larger scales gives the same form as Eqs. 11, but with 
scalar parameters which will have an explicit dependence on space, frequency and the values of the field in general[26]. 
The rate of spatial variation of these parameters at scales of the order of 50 pm is approximately zero, such that: 

' V < trj > | 10 6^ 3 ~ 

(12) 

, V < e w > 1 1 6 pm 3 ~ 

Note that there can be a frequency dependence of the current propagation which results from microscopic inhomogeneities 
of the electric parameters [6], from polarization phenomena [10] or from ionic diffusion [7]. This frequency dependence 
of the current will not disappear when considering larger scales. On the other hand, new frequency dependencies may 
appear, such as for example the transformation of a frequency-independent conductivity tensor at small scales (~ 1 pm) 
to a scalar conductivity at large scales (~ 50 pm) will be associated to a frequency dependence of this macroscopic 
conductivity (for details, see [7]). In agreement with this, measurements of the macroscopic conductivity demonstrated 
strong frequency dependence in different biological tissues [12]. 

E. The linear approximation in mean-field 

Still within the quasistatic approximation, we now consider the further simplification that the extracellular medium is 
linear. In a linear medium, the electric parameters are independent of the values of the fields (note that this linearity is 
different than that of Maxwell equations, which are always linear). In this case, the electric parameters only depend on 
space and time (or space and frequency). 

In this case, the system of equations (10) becomes linear at small scales. This linear approximation is easy to justify 
for the magnetic field, given the small currents involved (for example, An x 10~ 7 H/m in neocortex) and the gradient of 
p is almost zero (see [13]). In contrast, the linear approximation is less trivial in the case of the electric field because 
of the many nonlinearities involved. For example, several ionic conductances are strongly voltage-dependent (such as 
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the Na + /K + conductances involved in action potentials), which will make the electric parameters of membranes strongly 
dependent on the electric field. Nevertheless, the total volume of tissue occupied by membranes is small compared to 
other regions where the linear approximation is valid, so biological tissues can in general be considered as linear. Note 
that this linearity is evident for low frequencies (<10 Hz), but it is less evident for high frequencies (>100 Hz), where 
nonlinear phenomena such as action potentials can have a major contribution. 



F. The Kramers-Kronig relations under the linear approximation 



As discussed above in Section II B, the causality principle determines a supplementary constraint on the possible values 
of tensors < fi >, < e > and < cr e >. In the linear approximation, one can show that, in general (for isotropic media of 
class Cz), the linking equation between the electric displacement and electric field takes the following form: 



f 

Jo 



D(x,t) = E(x,t) + f(x,T)E(x,t-T)dT (13) 



In this case, one can show that the frequency dependence of electric parameters is not arbitrary but is linked by the 
Kramers-Kronig relations (see Section 82 in [14]): 

(14) 

tjJ y ' V/ n JO oi 2 -<i) 2 

where principal integrals (j-) are used. £«, is the absolute electric permittivity of vacuum and cr e is the static electric 
conductivity (<d = 0). Note that these relations can be seen as a direct and inverse transform. The Maxwell-Wagner 
time | | represents the characteristic time (or "inertia") for settling into a stationary regime, and can be strongly fre- 
quency dependent [10]. Interestingly, this ratio is mathematically analogous to the time-frequency uncertainty principle 
in Fourier transforms when the electric conductivity of the extracellular medium becomes very small at zero frequency 
(see Appendix A 5). Note that there exists no such relation for the spatial variations of electric parameters, which are 
specific to each medium. 

It is important to note that the Kramers-Kronig relations have a strong consequence on the plausibility of purely resistive 
media. If a medium is purely resistive, then both conductivity and permittivity are constant and independent of frequency. 
However, if one takes into account a very weak frequency dependence of conductivity ("quasi-resistive" media), such as 
~u)~°~() ~ f" witha << 1, then the Kramers-Kronig relations impose that we necessarily have fi^-Soo ~ f-^-°\ In such 
a case, the permittivity will be strongly frequency dependent, so will be the Maxwell-Wagner time r MW = |^|. Thus, a 
purely resistive extracellular medium is a singularity, and is not likely to be realistic model for complex biological media. 



III. APPLICATION OF THE QUASISTATIC MEAN-FIELD THEORY TO LINEAR MEDIA 



In this section, we apply the theory outlined above to media which are linear in the electromagnetic sense, and which 
are homogeneous and isotropic (at macroscopic scales). This is equivalent to assume that the macroscopic parameters 
< cr M >, < > and < D > are scalars independent of space, such that V(< cr u > +ito < e M >) x and V < D >« 0. This 
approximation is certainly valid for relatively large distances (greater than ~ 50 ymi). We also consider the system under 
the quasistatic approximation as defined above. 

In the following sections, we examine different limit cases. The first case corresponds to the standard model with 
dipolar sources and a resistive (or quasi-resistive) extracellular medium. A second case will consider the same model, but 
with additional monopolar sources. The third case will consider ionic diffusion (not present in the two first cases), which 
also implies monopolar sources. In each case, we will derive the expression to use for CSD analysis. 
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A. Dipole sources in resistive and quasi-resistive media 



We start with the standard model in which the electric conductivity of the extracellular medium is constant (in space 
and frequency), scale invariant, and isotropic. We also consider that diffusion is negligible. Under these hypotheses, we 
have V < y u >= (homogeneous and isotropic medium) and < y„ >= cr — est., and the general formalism (Eq 11) 
reduces to: 

v 2 < y„ > =o (15) 



We can also consider a slightly more realistic model of the extracellular medium by assuming that it is quasi-resistive 
(< cr m >« cr) instead of resistive (< <x w >= cr) (because the latter represents a singularity as outlined above in Section II F), 
isotropic and homogeneous for large scales (~ 50 fim; the medium is allowed to be non-homogeneous for smaller scales 
of the order of ~ 1 fim). Note that in general, the non-homogeneity of conductivity at smaller scales can induce a 
frequency dependence at larger scales (see [6, 7]). Thus, the hypothesis that the macroscopic conductivity is independent 
of frequency is equivalent to assume that there is no significant variations of impedance at microscopic scales. 

In such conditions, Eq. 1 1 again reduces to: 

V 2 < V l0 > = (16) 



In some formulations, the standard model does not consider the possibility of microscopic (~ 1 fim) monopolar 
sources [1, 3]. This is equivalent to hypothesize that, at every time, each portion of cell membrane has an equal number 
of positive and negative charges at opposite sides of the membrane, such that it is locally neutral. This hypothesis is also 
equivalent to state that the extracellular electric field is produced by dipoles (or more complex multipolar arrangements), 
and that the monopolar component of the field is negligible at scales of ~ 1 //m[27]. This implies that the attenuation 
of the extracellular potential follows a law which varies as 1/r 2 (or 1/r 5 , 1/r 4 ... for multipoles of higher order) when 
r — > oo, where r is the distance to the source. Thus, in the standard model, the electric displacement in frequency space 
(see Eqs. 4 and 5) is given by: 



< d:, > 



, < Eoj > + < Poj > + < 6o> > = < b u >< E w > + < d, 



(17) 



where 



V- < D* > = 



V- < P„ > 



, ^ n bound ^ 



^ Joound ^ 
Flo 



(18) 



V- < C a > = 

at large scales (~ 50 /mi). Taking the inverse Fourier transform, one obtains: 

V 2 < V > =0 



(19) 



According to this model, the inverse solution (CSD method) can be obtained assuming that the voltages measured at 
n different extracellular sites are solution of Laplace equation. According to the superposition principle, the extracellular 
potential can be considered as resulting from a sum of n macroscopic dipolar sources for sufficiently large n. Note that 
the value of n is determined by Shannon's sampling theorem, according to which the number of samples (number of 
electrodes «) must be twice larger than the larger spatial frequency of the field. To evaluate these n dipolar sources, one 
can simply apply the inversion of the matrix linking the n measured voltages with the n dipolar sources according to the 
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9 



"forward" solution of Laplace equation (see e.g., [3]). Note that this approach is different than the classic CSD method 
proposed by Mitzdorf [1], which is based on a Poisson type equation. 

The hypothesis of local neutrality in a homogeneous and isotropic extracellular medium implies that the frequency 
dependence of the measured signal is only due to the frequency dependence of the source (for example the effect of 
morphology - see [15], the exponential or bi-exponential nature of synaptic conductances, correlations in synaptic activity, 
action potentials, etc), because Laplace equation does not explicitly depend on frequency. Thus, according to the standard 
model, there is no filtering due to extracellular space and the power spectrum of the extracellular potential is identical to 
that of the current sources. 

Finally, it is important to note that, in a resistive extracellular medium, if we express the extracellular potential as 
a function of the dipole moments instead of the current sources, then the power spectral density (PSD) of the electric 
potential will necessarily have a supplementary frequency dependence of the form 1 /o> 2 compared to that of the current. 
This is due to the fact that the current is proportional to the temporal derivative of the the dipole moment p M (defined from 
the charge distribution). However, the situation is different if the medium is quasi-resistive. In this case, the Kramers- 
Kronig relations give e w ~ ^, and thus the ratio j± will have little frequency dependence. It follows that the PSD of the 
extracellular potential will have approximately the same frequency dependence as the current sources in a quasi-resistive 
medium. This is a striking difference between resistive and quasi-resistive media. As discussed above, the latter is a 
more realistic situation because any spatial variation of microscopic conductivity will necessarily induce a frequency 
dependence of the macroscopic conductivity (see [6]). 



In the previous section, we hypothesized that the extracellular medium is locally neutral at microscopic scales (~ 1 //m), 
and thus, that the sources of the electric potential are dipoles. We now relax this hypothesis, and allow significant electric 
monopoles to appear in addition to conductance variations, so that the field results from both monopolar and dipolar 
contributions. Electric monopoles could result from different physical sources, such as the ionic selectivity of synaptic 
ion channels (similar to a "Maxwell Daemon"), combined with the finite velocity of charge movement [16]. These factors 
should create some accumulation of charge in the vicinity of the synapse when synaptic conductances are activated[28]. 
Note that monopoles are transient by definition, and equivalently, one could consider that the conductance variations 
determine a non-stationary regime V- < j*> +— Kr — = (see Appendix A 2). In this transient regime, Kirchhoff's "point 
rule" does not apply (it is based on the law of current conservation V • j* = 0) and would apply only when the system 
reaches a stationary state. However, Kirchhoff's "loop rule" is always valid under the quasistatic approximation, because 
we have V x = 0, and consequently j> E ■ ds — 0, which is at the basis of the latter rule. 

Contrary to the assumptions of the dipole model, monopolar sources imply that integrating the charge density over a 
closed surface surrounding each source is non-zero. To include the contribution of monopolar current sources, we have 



where Q w the total charge contained in the source. Note that this relation shows that the monopolar component is linked 
to the current through a temporal derivative, which is a consequence of the charge conservation law. Consequently, the 
extracellular potential (which is here proportional to the charge) will not have the same power spectrum as the source, and 
will have an additional ~ -K component for a resistive extracellular medium. However, similarly to the case of dipolar 
sources in previous section, the situation is different for a quasi-resistive medium. The Kramers-Kronig relations imply 
~ i, and the ratio ^ will have very little frequency dependence and the PSD of the extracellular potential will be very 
similar to that of the sources. 



B. Monopolar sources in resistive and quasiresistive media 




(20) 
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If we consider the same conditions as for the standard model (resistive or quasi-resistive media), we obtain 

- v<di> = <p£r> 

v • £oo < g u > = < P iz e > + <pi " md > + < P h ° und > 

(21) 

V- < P m > = - < p h ° md > 

. v- <d> = -< P t cond > 

In such conditions, Eq. 1 1 becomes: 

V 2 < V w > =0 (22) 
for resistive and quasi-resistive extracellular media. Thus, in temporal space, we have the same equation for both cases: 

V 2 < V > =0 (23) 

However, if we take into account monopolar current sources and the law of charge conservation, then we have in 
general: 

cr„V 2 < Vuj > = i(L><p u> > (24) 

where cr^ = est. 

Thus, the model with monopolar current sources has a different structure than the dipole model in Section III A because 
< Po> >^ 0. Local neutrality in a homogeneous and isotropic extracellular medium implies an identical frequency depen- 
dence of the current source I w = -iu> < p w >) and the extracellular potential. Like the standard model, there is no "filter" 
in this case. There is a notable difference, however. The law of attenuation with distance varies here in l/r instead of 1/r 2 
for r — > oo. If the number of electrodes is large enough to respect Shannon's sampling theorem, then the current source 
densities can be simply evaluated by approximating the Laplace equation using finite difference methods, as well as the 
knowledge of the "forward" solutions of this equation (see [3]). We will see in the next section that these conclusions are 
different if ionic diffusion is taken into account. 

C. Models with ionic diffusion 

While the influence of ionic diffusion was neglected in the previous sections, we now consider this case more explicitly 
without any other hypothesis about the medium. If a selective ion channel opens, the flow of ions may induce accumulation 
of charges in the region adjacent to the channel if ions diffuse faster than the time needed for passing through the channel 
(which will generally be the case). The electric field resulting from conductance variations is not selective on the type 
of ion, such that the positive ions are attracted and negative ions are repulsed if the field is negative (and vice-versa for a 
positive field). This is contrary to the flow inside the channel because it is selective to only a subset of ionic species. The 
combination of these factors makes it unavoidable that there will be charge accumulation around open ion channels. In 
the standard model, this charge accumulation is considered as negligible. 

We now evaluate the consequences of this phenomenon on the frequency dependence of the field produced by ionic 
conductances in the subthreshold regime. If we consider a homogeneous extracellular medium with constant electric 
parameters (independent of frequency at large scales, ~ 50 pm), then we have: 

< <r£ > \ M = a- 

< Soj > \ M = s (25) 



<D>\ M = D±0 
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where the parameters cr, s and D do not depend on space. 

According to those hypotheses, variations of ionic concentrations appear in the vicinity of the open ion channels, 
and these variations are opposite to the current produced by the electric field resulting from conductance variations. It 
thus appears that the conditions of current propagation at microscopic scales (~ 1 /mi) cannot fulfill the condition of 
homogeneous ion concentration which is at the basis of Ohm's differential law (see Appendix A 3). In such conditions, 
the electric parameters of the extracellular medium have the following form at microscopic scales (~ 1 /mi): 

< cr* > \ m (x) = &„,{x) 

<So>> \m(x) = Sm(x) (26) 

< D > \ m (jt) = D m + 

with 

< 7a> > \m =<<y,>> \ m > \m (27) 
This last equation is necessary to keep the consistency between microscopic (~ 1 pm) and macroscopic (~ 50 /mi) scales. 
According to this model, the current density (at microscopic scales, ~ 1 /mi) is given by: 

<J> \m = -<?mV < V > \ m + D m V < p > | m (28) 

where we have (see Eq. All) 

o- m = A q T c n v (x, t) (29) 

This expression can be deduced by separating the domain into sufficiently small elements such that ion density can be 
considered as spatially homogeneous, and sufficiently large for Ohm's law to apply. 

According to Boltzmann distribution (see Appendix A 1), we have 

2 2 

V < p > u*, t) = «„V < V > \JX t) = ^-i^-o- m V < V > \ m (30) 

By taking into account Eqs. 28 and 29, we obtain 

-> - kA„T r T 

< j > U = [Dm - — d v <p > | m =</? > L V <p > L (31) 

< q > L 

where < /3 > is an effective diffusion coefficient. Note that the value of < ft > is smaller than the mean diffusion 
coefficient because J T ^ 2 must be positive. The value of < [3 > also depends on the values of ionic concentrations because 
several parameters in Eq. 31 are concentration-dependent and is proportional to temperature because the ionic diffusion 
coefficient is itself proportional to temperature (see for example the Einstein relation for diffusion). 

Applying the differential law of charge conservation, we get 

</j>| m V 2 <p>|,„^- g< ^ >1 '" (32) 

ot 

Thus, the charge density produced in the vicinity of the ion channel is solution of a parabolic differential equation similar 
to the diffusion equation [29]. 

It follows that the charge density obeys: 



V 2 < poj > |,„ = -i <Pco> \m (33) 

< P > \m 



to 



At microscopic scales (~ 1 /mi), we obtain (see Eq. 11): 

vr2 „ , V(<y<o>\m) _ I7 , 10) <D>\„, . 

V 2 < > \ m + ■ V < Voi > U = ■ — — < p w > \ m ~iuxV > \ m (34) 

< 7o> > \m <7o,> \m <P> \m 
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Here, the proportionality between < p w > \ m and < V w > \ m can be deduced from the linear (first-order) approximation of 
Eq. A3 (see Appendix A 1). The second-order approximation would give a cubic term in < V w > | m . 

Applying the consistency equation between scales by assuming the statistical independence of the parameters leads to 
the following equality: 

1 

< Ju > \M <P> 

with 

N 

1 - 1 U 



9 lUJ 

V z < Voj > \ M = 1 — — < Pco > \m (35) 

\M 



<7oj>\m N p l <yi>\„ 



</3>\m N ^ <jP>\, 



<D>\ M = j; 2 < DJ > 



1 N D' | 

N 2 <fi/>|.„ (36) 



where N is the ratio between the reference volumes at macroscopic and microscopic scales (note that to simplify the 
formalism, we have approximated the macroscopic mean by a discrete summation over microscopic means). The second 
term of the lefthand side of Eq. 34 becomes zero at macroscopic scales (see consistency equation Eq. 27). Note that the 
means over parameters y w and /3,„ are harmonic means, while the means over matter fields are arithmetic means. 

Finally, by applying the inverse Fourier transform, we obtain (for Maxwell- Wagner times much smaller than unity): 

crV 2 <V>| M = J—?L <p> \ M (37) 

<P> \hidt 

Thus, the CSD method in the presence of ionic diffusion takes a form which is very close to the Mitzdorf model [1], 
because we have one source term. However, there are two notable differences: first, the frequency dependence of charge 
density implies that the extracellular medium will be frequency dependent according to an impedance which varies as 
1 / V(w) (see Appendix A 4). Second, the extracellular potential attenuates with distance according to a Yukawa potential 
- — instead of 4-, as in the standard model. In this case, we have 



<Vo>> | m (r)| = \<Vcj> \ m (R) | (38) 



Here, the extracellular potential is proportional to the charge density (see Eqs. A3, A14, and A15 in Appendix A 4) 
under the linear approximation and for a spherical source. It is interesting to note that the exponential term increases 
with frequency such that the extracellular medium favors the propagation of low frequencies (low -pass filter), as shown in 
Fig. 1). This type of attenuation law in Fourier space is also consistent with an exponentially-decaying impedance. If the 
frequency spectrum is narrow, one can replace u> by its maximal value k — | ^ ''^" ) |, which leads to an attenuation law 

for the potential as e ' r f °' . We can thus write that the electric field is approximately equal to < E > \ m = -V < V > |,„ = 
e^ k ^ r " > (kr+\)-Kf. In a resistive medium, this leads to an electric resistivity given by < p e >„,= <o _ f ^ (r) = ^fe ~ k(r ~ rtt) (kr+l) 
(note that the current / = \j\ 4nr 2 is conserved in a resistive medium, and the field is given by < E > \ m -< j > \ m j < 
o~( r ) > D- Note that this particular distance profile of the potential was calculated by assuming that the medium is 
homogeneous (see Appendix A 4), which makes it applicable only at short distances from the membrane (of the order of 
10 to 50 nm). 

[ Figure 1 here ] 

Finally, it is important to note that this frequency dependence cannot be removed because it is an effect of the feedback 
caused by ionic diffusion when ion channels open, and this is inherent to biological tissue. Because the PSD of the 
extracellular voltage is of the form ~ ^ I(u>) (see Eq. A15 in Appendix A4), one can view the effect of ionic diffusion as 
a "1/f filter", as found previously [7]. 
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We now compare the different cases examined here. From the point of view of the differential equations involved, in 
the "standard" model based on dipoles, as well as with monopoles, the extracellular potential is solution of the Laplace 
equation, which is elliptic. In the third model with ionic diffusion, the extracellular potential is solution of a Poisson 
type equation where the source term is proportional to the time derivative of the voltage (under the linear approximation), 
which gives a parabolic equation. As outlined above, the diffusion model is closer to the monopole model (as diffusion 
can have monopolar effects) in a resistive extracellular medium, but leads to a fundamentally different mathematical form. 
The physical reason for this difference is that the ionic diffusion at the interface ion channel/medium increases the inertia 
of the system as a function of frequency. 

At the point of view of the CSD analysis method, different algorithms must be used according to which model of the 
extracellular medium is assumed. In the two first cases, one must use a "forward" solution because Laplace equation is 
non invertible. In this case, it is necessary to explicitly include the distance dependence of the extracellular potential , 
which varies as 1/r 2 for dipoles (for distances sufficiently large compared to the size of the dipole) and 1/r for the model 
based on monopoles. Note that if the distance to the sources is not large enough (compared to the typical size of the 
sources), or if the dipolar moments are very large compared to monopolar moments, then the attenuation will be closer to 
a linear combination of 1/r and 1/r 2 . 

In the diffusive model, however, the approach is totally different because of the parabolic nature of the equations. In this 
case, it is enough to apply the Laplace operator to recover the sources. Two strategies are possible. First, one could simply 
apply Laplace operator on the extracellular voltage to yield estimates of the current source densities. Second, one could 
use a "forward" model and consider an attenuation law following a Yukawa potential e M)r and apply the same procedure 
as for the other models. 

Perhaps the most interesting aspect is that the three different models considered here have a different spectral signature. 
In the dipole model (or monopole model in a quasi-resistive medium), the PSD of the extracellular potential is identical 
to that of the sources. The resistive monopole and the resistive dipole model exert a filtering effect of l// 2 type, whereas 
the diffusive model is equivalent to a 1 // filter. Thus, the frequency characteristics of the signal can serve as a criterion 
to determine the most appropriate model. For example, if the PSD of the extracellular voltage has l/f structure, this 
automatically discards a pure monopole model, as well as dipole models in resistive or quasiresistive media, and would 
suggest diffusive type models. 

E. Synthesis and applications to experimental data 

In this section, we synthesize the theoretical developments provided here, and suggest a guide of how to apply them to 
experimental data. The generalization of the CSD method for different cases of current sources and type of extracellular 
medium, is summarized in Table 1 . The table considers monopolar and dipolar current sources, as well as different types 
of resistive and non-resistive media. 

To perform a CSD analysis by allowing non-resistive properties of the extracellular medium, we suggest the following 
procedure. 

1. Estimate the type of extracellular medium from the power spectral structure of CSD signals. As detailed above 
(Sections III A, IIIB, IIIC, the type of medium (resistive, quasi-resistive, diffusive, etc) and type of current sources 
(monopolar, dipolar, etc) can be inferred from the power spectral structure of the extracellular potential. This 
analysis should be done on non-filtered data to set constraints on the possible combinations of sources type of 
medium. 

2. Identify the correct CSD expression compatible with the type of source/medium inferred from power spectra. Ta- 
ble 1 summarizes the different cases considered here. The expression identified is then used to calculate the current 
sources from the extracellular potential recordings. 
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TABLE I: Different generalizations of the CSD method. The table shows the mean-field equations for different types of media, and 
for monopolar or dipolar sources. The Mitzdorf model is shown apart, because it does not correspond to any of these mean-field 
scenarios. Abbreviations: res. => resistive homogeneous medium, quasi-res. => quasi-resistive homogeneous medium, res. + dif. => 
resistive homogeneous medium + ionic diffusion, gen. => general (non-homogeneous, with spatial and frequency-dependent variations 
of electric parameters), 1-pole => monopole, 2-pole => dipole. Note that the frequency-dependence of the permittivity and conductivity 
are not independent but are linked by the Kramers-Rronig relations. The quantity ^ is the ratio between the Fourier transform of the 
extracellular potential V u and the Fourier transform of each point current-source which produce the field (asymptotic solution, far 
from the sources). The function f(co) in "res+dif" determines a Yukawa type potential (see Fig. 1). 

IV. DISCUSSION 

In this paper, we have formulated a series of generalizations of the CSD analysis method applicable to extracellular 
recordings in brain tissue. This generalization is based on a general theory that we derived and which aims at linking the 
extracellular potential with current source densities in the tissue. We have considered a mean-field version of Maxwell 
equations by considering the different fields as averages over some reference volume. By varying the size of this volume, 
one can apply the same theory to different scales. At microscopic scales (~ 1 |/m and smaller), the theory must use 
the microscopic values for electric parameters (for example, the very different resistivities of fluids or membranes). For 
mesoscopic or macroscopic scales (~ 50 //m and larger), the theory can directly include the "macroscopic" measurements 
of conductivity and permittivity, as well as their possible frequency dependence if needed. Note that this mean-field 
approach takes into account the physical and biological properties of the sources, and thus is more general than previous 
approaches [7] which did not consider source densities. 

We have examined different limit cases, such as a purely resistive extracellular medium with current sources consisting 
exclusively of dipoles, in which case the theory recovers the standard model. In this standard model, the mean-field theory 
shows that the electric potential must be solution of Laplace equation, such that the "classic" CSD approach of Mitzdorf [ 1 ] 
does not apply. To inverse the CSD in this model, one must apply the forward solutions of Laplace equation because the 
associated operator is non-invertible (see [3]). In resistive media, the extracellular potential must have an additional 
frequency dependence of 1 // 2 relatively to that of the current. Interestingly, we found that Laplace equation remains 
valid for extracellular media which are quasi-resistive (where the electric parameters weakly depend on frequency). In 
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this case, the frequency dependence of the extracellular potential is similar to that of the current. A weak frequency 
dependence was indeed found in some experimental measurements of resistivity [17, 18], while other experiments [12] 
displayed a much more pronounced frequency dependence. With respect to the attenuation with distance, the standard 
model predicts an attenuation law as 1 /r 2 , for both resistive or quasi-resistive media. 

We also examined the case of monopolar sources. If such monopolar sources are present in addition to dipolar sources, 
within resistive or quasi-resistive media, then the CSD equation takes a slightly different form predicting that the potential 
will attenuate asymptotically with the inverse of distance (l/r), while the standard dipole model predicts a square depen- 
dence (1 /r 2 ). With monopolar sources, the potential in the extracellular medium is also solution of Laplace equation, and 
thus the inverse algorithm of the CSD method does not apply identically as for dipoles. To find the inverse CSD, one 
proceeds similarly as the standard model by using the "forward" solution of Laplace equation. However, in this case the 
sources must be considered as a linear combination of terms varying as 1 jr (monopoles) and l/r 2 (dipoles) in this forward 
solution. 

As a third model, we examined the case of ionic diffusion within resistive or quasi-resistive media. In this case, the 
CSD takes a form very close to the "monopolar" CSD discussed above, but we found that charge density is frequency 
dependent according to a Warburg impedance in 1 / V(w) (see Appendix A 4 and Section III C). This result is in agreement 
with a previous modeling study of extracellular potentials in the presence of ionic diffusion [7]. Another consequence is 
that, for spherical symmetry, the attenuation with distance follows a Yukawa potential e , which decays faster than the 
different laws considered above for large enough frequency (see Fig. 1). This particular form is responsible for a low-pass 
filtering of the extracellular medium. Note that this form is obtained for spherical symmetry, but other forms may be 
obtained in different geometries. 

It is important to note that the CSD theory was originally designed without specific hypotheses about the nature of cur- 
rent sources [2, 4, 5], other versions of the CSD theory clearly assumed that current sources are dipoles [1, 3], Assuming 
dipolar sources is equivalent to assume that we have stationary current conditions at all scales. However, we show here 
(Appendix III B) that at small scales (synapses), such a stationary current condition is not necessarily met. A first possible 
source of monopolar effects is the inertia of charge movement along membranes together with ion-channel selectivity. 
Following the opening of ion channels, the flow of ions will entrain a re-equilibration of the charges adsorbed on both 
sides of the membrane. While this process is usually considered as instantaneous, together with neglecting ion-channel 
selectivity, these processes may have important consequences. Indeed, if one takes into account the fact that charges 
do not move instantaneously and ion-channel selectivity, this will necessarily create transient charge accumulation and 
monopoles. A similar effect will occur through ionic diffusion and electric field, at the interface between the ion channel 
and the extracellular medium, because ions diffuse faster than their mean passage time through the channel, which will 
also create charge accumulation and monopolar effects. Note that when this mechanism produces an external electric 
field which will contribute to the extracellular field (in addition to transmembrane currents). These effects will contribute 
to transient monopoles, during which Kirchhoff's node law will not apply. The fact that ions move considerably slower 
than electrons in a metal conductor will also participate to deviations from Kirchhoff laws. Whether this transient time 
is significant, and whether the system could be continuously "outside of equilibrium" due to sustained synaptic activity, 
should be investigated by future work. 

In the diffusion model, one can directly use the Laplace operator to inverse the CSD, contrary to the other models. 
Taking into account ionic diffusion requires to revise the "forward" approach, because the attenuation law does not follow 
a l/r 2 law, but rather a Yukawa-type law while the extracellular medium is associated to a Warburg-type impedance. 
In a previous study, we showed that indeed a Warburg type impedance could account for the transfer function between 
intracellular and extracellular potentials [19] (for frequencies comprised between 3 and 300 Hz). It is also consistent with 
measurements of conductivity and permittivity [12] (but see [17]). Note that the linear approximation in the diffusion 
model is not valid for high values of the potential (larger than ~50 mV; see A 4), so this model applies well to subthreshold 
activity, but may need to be revised for action potentials. Similarly, corrections to the CSD given by the "forward" 
approach (see for example [21]) may also need to be reformulated for non-resistive media. 

Thus, with respect to the paradox of the CSD method, as described in the introduction, our study suggests that it is 
naturally solved by taking into account ionic diffusion. This introduces an additional source term in the general equation 
for the electric potential (see Eq. 11). This additional term gives a Poisson type equation for the potential (instead of 
Laplace equation), similar to the classic CSD approach. Contrary to the cases with resistive and quasi-resistive media, the 
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classic algorithm of CSD inversion given by Mitzdorf [1] is applicable here. Thus, the results obtained with the classic 
CSD analysis are perfectly consistent with ionic diffusion because diffusion gives a source term which is very close to 
the phenomenological model of current source density introduced by Pitts and Mitzdorf [1,2], but in a manner consistent 
with Maxwell-Gauss law. So, we conclude that the usual approach for CSD inversion, although paradoxical, should 
nevertheless give results equivalent to a model with ionic diffusion and consistent with Maxwell-Gauss law. 

Finally, the few limit cases considered here are by no means exhaustive. For example, we neglected the Maxwell- 
Wagner time of the extracellular medium and the microscopic variations of impedance. The theory outlined here is 
general enough to include these effects if needed, which is another way to solve the paradox. For instance, considering 
phenomena such as "reactive" extracellular media, which react to the electric field (for example through polarization of 
cell membranes), can be done by taking into account the Maxwell-Wagner time of the medium (see details in [7, 10]). 
According to Gabriel et al. [12], the macroscopic electric permittivity becomes larger while macroscopic conductivity 
becomes smaller for smaller frequencies, when the electric field is imposed according to a well-defined direction. In 
these measurements, wr MW << 1 for frequencies larger than 10 Hz, but wt mw may be considerably larger for lower 
frequencies [7], where electric polarization may play an important role. The second term in the lefthand side of Eq. 11 
would then not be negligible anymore. Because this term can be considered as an additional source term (see [6]), similar 
to the case of diffusion, this also solves the paradox described in the introduction. 

In conclusion, we have provided here a generalized CSD approach valid for more realistic properties of the extracellular 
medium, taking into account ionic diffusion or polarization effects, usually neglected in the standard CSD analysis [1, 3]. 
We found that including such effects may have deep consequences on the expression to be used for estimating current 
sources, and thus may also have consequences on the values of current sources estimated from experimental recordings. 
For example, the potential due to monopolar sources will decay slower than for dipoles, which will necessarily affect the 
recorded potential at the electrode. Similarly, considering "reactive" aspects of the extracellular medium by including a 
significant Maxwell- Wagner time leads to a different CSD expression, close to the form derived for ionic diffusion. Future 
work should apply these expressions to extracellular recordings in brain tissue, with the aim of identifying which of these 
phenomena are most consistent with experimental data. 



In this appendix, we consider ionic diffusion at the interface between ion channels and the extracellular medium, as well 
as at the interface with the cytoplasm. We use the quasistatic approximation (in the thermodynamic sense), which implies 
that the net charge density must be solution of a parabolic partial differential equation, as for pure diffusion phenomena. 
We will next consider system in spherical symmetry, in which case the impedance is equivalent to a Warburg impedance. 



Appendix A: Appendix 



1. Impedance for systems with ionic diffusion 



a. Ionic diffusion under the quasistatic approximation in the thermodynamic sense 



The current density at microscopic scales obeys the equation: 



< j > |,„ = -cr m V < V > |„, + D m V < p > I 



\m 



(Al) 



with (see Eq. All) 



= A q T c n D (x, t) 



(A2) 



Let us assume that the system is in a quasi-static case in the sense of thermodynamics. As shown by application of 
Maxwell distribution of velocity distribution and the principle of detailed balance [20], we can deduce the Boltzmann 
distribution for a field which varies infinitely slow. This approximation is valid here because the drift velocity of ions 
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under an electric field is much lower than the absolute velocity of ions (which is of the order of sound velocity). Within 
this quasistatic approximation, we can apply the Boltzmann distribution to obtain the number of ions per unit volume as 
a function of time and space: 

_, <,>i„,<y>[„,(?. Q < ? >l,„<y>l„,(?. t) 

n D (x, t) - n™ [e + a + e a ] 

when we assume that V(oo) = at infinite distance, where njj° is the number of ions per unit volume at an infinite distance 
from the source ("far distance") and < q > is the mean absolute charge, k = 1.3806503 x 10~ 23 J/ °K is the Boltzmann 
constant and T is the temperature in degrees Kelvin. It follows that the net charge density is related to the value of the 
electric potential according to: 

<p>UKf)=rC<9> L [e + * e n ] (A3) 

where < p > |,„(n, t) is the average net charge density. Note that the - sign in the second righthand term comes from 
the sign of the charge. Also note that this relation implies that the net charge density is zero at an infinite distance, and is 
linked to the electric potential by a nonlinear relation[30]. 

Applying the operator V on the net charge density gives: 

v < P > | m (f, t) = (<q l > J m)2 «„v < v > \ m (x, t) = {<q> }"f cr m V <V> \ m (A4) 

kl AqT c kl 

Note that no such relation would be possible outside of the quasistatic approximation in the thermodynamic sense. 



2. Non-stationary aspect of the electric field produced by membrane conductance variations 

In this appendix, we show that if the transmembrane current is non-zero, this necessarily implies a charge variation at 
the interior of the compartment according to a non-stationary regime 

J dt 

We define as "interior" the domain delimited by the inner surface of the cell membrane as indicated in Fig. 2. We will 
show that for different values of stationary current, the net charge inside a cable compartment is different from zero, and 
depends on the value of the current. Thus, when the current is variable, the system is necessarily in a non-stationary 
regime. 

To demonstrate this, suppose that we have a current density /which is time independent and stationary, satisfying 

V-/=0. 

The amount of charge situated in the interior of the compartment strictly inside the compartment can be calculated from 
the integral of the electric field across a closed surface: 

Qi„t = e£^ ■ hdS — JJ" s c ■ h dS + JJ" e m E ■ h dS (A5) 

3D channel+axial membrane 

where the integral is made over a Gauss surface which goes through the middle of the membrane thickness, as indicated 
in Fig. 2. The surface also avoids ion channels (by surrounding them below their inner side), e = e c and e = s m are the 
electric permittivity of the cytoplasm and of the membrane, respectively. 

Taking into account Ohm's differential law and stationary current condition, one obtains: 

/ = Jj j-hdS = jjj V-fdv= jj cr c £ ■ h dS = (A6) 

channel+axial channel+axial 
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where cr c is the electric conductivity of the cytoplasm. Because the capacitive impedance is infinite for / = 0, the surface 
integral over the membrane is zero, and we obtain: 




(A7) 



membrane 



because s c e! = when electric parameters are independent of space. 

The cylindric symmetry together with the isopotentiality of the surface of the membrane compartment imply that this 
integral is non zero and depends on the value of the transmembrane current /,„ . If we have two different transmembrane 
currents, then the values of the electric field inside the membrane are necessarily different and the values of Q lnt are 
different too. In particular, if the current /,„ is zero, this integral equals the negative charge that there can be inside the 
compartment at rest. 

It follows that, if we have a variable transmembrane current, the interior charge varies to satisfy ^ ' j - ~% ^ because 
the membrane time constant t„, is not negligible. This shows that the variation of current sources caused by membrane 
conductance variations have the properties of a non-stationary regime. Note that this non-stationarity requires to take into 
account the volume of the membrane compartment, and would not be present for point processes. 



3. Ohm's law and frequency dependence 

In this appendix, we show that Ohm's law implies that the ratio between potential and current does not depend on 
frequency for frequencies smaller than 10 4 Hz. The aim of the appendix is to show the physical bases of the model shown 
in Section IIIC. 

Let us apply an electric field E, independent of time and space, at time t — and during a time interval AT, to a 
homogeneous aqueous solution (similar to salted water) at thermodynamic equilibrium. The field will accelerate every 
ion according to the following velocity law: 

%) = m + —E {ff - td (A8) 
m q 

where f,- is the initial ion collision time and tf its final collision time. If the field is applied during a time interval AT much 
longer than the typical collision time between two molecules (or ions), then the mean velocity of ions is given by: 

< v(t) > A ;r=< v{ti) >&t +^-3 < tf - t, > A7 -=< v(ti) > AT +-^r£r* (A9) 
m\ m\ 

where the mean collision time of ion k, qi the charge of ion k and the mass of ion k. Note that the mean collision 

time for ions such as Na + , K + , Cl~, Ca ++ in sea water are approximately equal to 10~ 12 10~ 14 s for temperatures 

between 250 °K and 350 °K [20]. 

If the system is at thermodynamic equilibrium, we have < i>(f,) >,= 0. It follows that the time average of the current 
density is given by: 

2 

where n k D is the number of ions of type k per unit volume. Because the quantities %, t* and «* are statistically independent, 
we can write: 

cr e - A q T c n D (All) 

where A q -< ^ >k, t c —< t£ >k and n„ - N < n\ >t are average values of the quantities t c and n v for each type of 
ion. 

Thus, this expression shows that the ratio between the time average of the current and the electric field does not depend 
on time when the time interval AT is much longer than the mean collision time between two molecules of the solution, 
because the number of collisions of each ion is very large. 
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Because the mean collision time is of the order of 1CT 12 - 1CT 14 s for a temperature included between 250 °K and 
350 °K (see details in [20]), we can write that we have a very good approximation of a variable electric field by a 
piecewise constant function (staircase), with time intervals are much longer than the mean collision times. Then, the time 
average of the current density will be proportional to the applied field. We can therefore write that the ratio between mean 
current density and applied electric field does not significantly depend on the frequency of electric field - nor of time - 
for frequencies smaller than 10 4 Hz [31]. This approximation is called differential Ohm's law. 



4. Ionic diffusion in spherical symmetry 

To simplify the computation of Eq. 33, we consider a model with a spherical source, surrounded by a spherically- 
symmetric (or isotropic) extracellular medium with boundary condition and isopotential over the surface of the source. In 
this case, we have the following equality in Fourier space: 

d 2 <p bJ > | m 2d<p 0J > L. u 

d* + 7 dr +l Tp^ t <^>l<»=° < A12 > 

The general solution of this equation is: 



<Pu> \ m (r) = A(cj) + B(cS) (A13) 

r r 

where r is the distance between the geometric center of the source. Because we must have no charge accumulation at 
infinite, we must set A = 0, which gives: 

< Po> > UW =< Poj > L(R) (A14) 



Taking into account Eq. 31, the current density is given by: 



d < Do, > \ m 1 I CO 

< jo > \m(r) =< /3 > U =< B > \ m ( - + J-i— — - ) < Poj > U(r) (A15) 

or r \/ <B> \ m 

By developing the net charge density (see Eq. A3) in Taylor series, we have at first order: 

<q>\ m <v>\„,g n «,>\ m <v>i,„(!. i) (< q > m ) 2 n°° 
<p> | ffl (f, = < < q > \m [e + * e f ] * 2 V H ™ v - <V> UZ t) (A16) 



At physiological temperature (310 °K), we can write that the precision of this linear approximation is larger than 90 % 
if the potential is smaller than 15 mV (these estimates are obtained by replacing Boltzmann constant by the values of the 
mean charge, ss 2 x 10~ 19 C). However, the precision drops to about 8 % for 100 mV, in which case one must consider up 
to the third-order term (~< V > |^,) in the Taylor expansion. This is the case for action potentials, and thus the propagation 
of the field will become more complex for spikes. 

It follows that the impedance between infinite distance and a given point P at distance r from the source is given by: 

kT 1 
ZJr) = 2- , o fl , co j (A17) 



Thus, the impedance will tend to a Warburg impedance for large distances from the sources, and for high frequencies. 
Moreover, if the curvature radius is very large, or for planar membranes, one can set R — oo, which gives: 
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where C ~ (1 - i) and thus, the phase of the impedance becomes independent of frequency. 

One sees that the value of the parameter < /3 > \,„ and the curvature radius determine the magnitude of the phase. 
They also determine the distance at which the impedance becomes equivalent to a Warburg impedance. For example, the 
estimates of impedance in rat cortex made previously [19] indicate that the impedance can be well approximated by a 
Warburg impedance for frequencies above 3 Hz. 



5. Kramer-Kronig relations 



In this appendix, we show that the Maxwell-Wagner time r MW = ~ {\ when the conductivity is very low at zero 
frequency, based on Kramers -Kronig relations (see Eq. 14). 

These relations show that if crj does not depend on frequency, then e w will also be frequency independent. These 
relations also show that if e u - Sco = koT^ \ then we have (see Eq. 14) 



rj - o-q = — — t " — da > 

71 Jo 



o co' z - co z 



By setting x = — , we obtain 



'to "0 



—2 f °° 

— T —r~, — s dx \ kto h — k\kco b 

L n Jo x l - b (x 2 - 1) J 

where we have by definition 

r° i r r 1_kl i f°° i i 

k\ - T — r~i — n dx — lim I ——. — z dx + I — -—. — z c/x' 

Jo ^(^-l) kHoLJ x 1 -*^ 2 -!) J 1+w x' 1 - b (x' 2 - 1) J 

By replacing x = l/x' in the second integral in the righthand term, we get 

pi-kl i 
h = lim [x^'^'+x^^]— rfx 



This relation shows that we have s " B ™ e = -J— = — (we have set k = j-) such that a small variation of <x w relative to 
frequency (for low frequencies) will entrain a strong variation of e w relative to frequency, and thus a very large value of 
t mw for low frequencies. If cr^ ~ and if e w is much larger than that of vacuum, then r MW is given by: 

_ £(o * /at n\ 

where /c is a constant which depends on the exponent b (see Fig. 3). This approximation corresponds well to the experi- 
mental measurements of Gabriel et cd (see [12]). 

Finally, if the conductivity varies very slowly with respect to frequency (small value of a), then the permittivity will be 
proportional to l/co, and will therefore vary very steeply at low frequencies. 
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Figure legends 



FIG. 1: (color online) Attenuation profile of the extracellular potential as a function of distance. The profile of the potential with 
distance is shown for two positive values of < p >: one value comparable to the diffusion coefficients of k + ,Na + , CI" (left), and another 
value 100 times smaller (right). When < /3 > is comparable to the diffusion coefficients, the attenuation according to Yukawa potential 
is similar to Coulomb's potential and we have a Warburg impedance. For smaller values of < ft >, the Yukawa potential determines a 
significant additional low-pass filter. The attenuation law is very steep for frequencies larger than 100 Hz when < /} > \ m = 10~" m 2 /s. 
The different curves indicated are the Coulomb's attenuation law as 1/r (monopolies; solid line) and as 1/r 2 (dipoles; dotted line; thick 
gray lines correspond to attenuation according to a Yukawa potential; thin lines correspond to Yukawa attenuation combined with a 
Warburg impedance, which is the most complete case taking into account ionic diffusion effects. In each case, different frequencies are 
compared (1, 10, 100 and 1000 Hz) and are shown by different colors and dashed lines, as indicated. The current source has a radius of 
5 run. 

FIG. 2: (color online) Scheme of an isopotential membrane compartment with ion channels. The blue circles indicate ion channels in 
the membrane. The dotted line indicates a Gauss surface delimiting the interior of the compartment. The values of electric permittivity 
e are equal to e„, in the membrane and e t . in the cytoplasm. The electric conductivity a is equal to cr c in the open channels, and is 
assumed to be zero in the membrane. 

FIG. 3: Representation of k as a function of the exponent b. The values of k were calculated from expression _ 2 ,„ — — — , where 

IT Tn 1-J,. 5 
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